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1 Introduction 

The correspondence between the two descriptions of 2D quantum gravity given by the matrix model (MM) and the 
Liouville gravity (LG) has been widely investigated for more than thirty years [1,2]. One of the key steps was the 
matching of the MM critical exponents with the gravitational dimension of the LG correlators [3-5]. Technically 
more elaborate, the matching between coiTclators was first performed in [6] for the sphere one and two points 
functions. This correspondence involves the so-called resonance transformation, a finite renormalization of the 
couplings arising due to the ambiguity lying in the presence of contact terms. Since these early steps, various matrix 
models have been introduced (one and two hermitian matrix models, ADE, matrix chains, 0(?i) loop gas model, 
q-Potts,...) to describe a discrete quantum gravity with different kind of matter fields. We concentrate here on 
the simplest one, the one hermitian matrix model which continuum description is provided by a minimal Liouville 
gravity, i.e. a quantum gravity which matter sector is given by a minimal model. In this context, the full resonance 
transformation on the sphere has been conjectured in [7], and checked up to the fifth order [7, 8] (see also [9] for the 
case of the gravitational scaling Lee- Yang model). Agreement between the two descriptions has also been verified 
on the disc with the trivial boundary conditions (BC) in [10]. 

Recently, the correspondence between the one matrix model and minimal LG has been extended to worldsheets 
with arbitrary boundaries. More precisely, matrix con^elators describing a disc with any boundaiy conditions have 
been constructed by two of the authors in [11]. It was later shown [12] that this construction relies on a linear 
relation of decomposition satisfied by the FZZT branes [13, 14] of LG found in [15]. In [16], we constructed the 
matrix boundary operators introduced between two arbitrary boundary conditions. This construction involves the 
first order of a pure boundary resonance transformation that takes the form of a polynomial in the matrix. ' 

The purpose of this paper is to extend the previous constructions by considering matrix coiTclators describing 
a disc with an insertion of one bulk and one boundary operator In the search for agreement between LG and MM 
correlators, we are led to distinguish two different cases. When the boundary operator Kac index does not belongs 
to the fusion rules of two copies of the bulk operator index, the matter part of the LG con^elator has to vanish. This 
case involves a mixed bulk-boundary resonance transformation that will be determined at the first order by imposing 
the cancellation of MM correlators combinations. On the other hand, if the LG correlator is non-vanishing, no 
mixed bulk-boundary resonance arises for dimensional reasons. Thus, we only have to consider the pure boundary 
resonance already determined in [16]. We explicitly checked the agreement between MM and LG in a few non- trivial 
cases, confirming our general approach to the resonance transformation. In this process, we derived an alternative 
expression for the polynomials encoding the resonance transformation. This expression proved to be convenient 
when the Kac index of the boundary operators are small. 

In the first section, we concentrate on the LG side and speciaUze the expression of the Liouville bulk-boundary 
coiTclator [17] to the case of a coupling with a degenerate matter operator. In section two, we explain how to 
construct the bulk-boundary matrix conelators and derive their expression in the continuum limit. The section three 
is the core of the paper, starting with general considerations on the resonance transformation and deriving the relation 
between MM and LG correlators. Then we show the agreement between the two approaches for both vanishing and 
non- vanishing LG correlators. Most of the technical details are assigned to the appendix. 

2 On the Liouville gravity side 

2.1 Bulk-boundary correlator in Liouville gravity 

The Liouville gravity is the description of 2D quantum gravity with critical matter in the conformal gauge [18]. The 
total action obeys the conformal symmetry and consists of three components, the Liouville, matter and ghost parts 



'strictly speaking, this resonance transformation also includes the bulk LG cosmological constant linked to the dressed bulk identity 
operator. In "pure boundary", we understand here that no other bulk operators are involved except for this cosmological term of the critical 
action. All these statements will be clarified in the section three below. 



The bulk Liouville action Sl is given as 



5l = ^ / ^(fz [g^'^dacpdbqy + QR(I) + ^irfie^^'^ 



(2.1) 



where fi is the bulk cosmological constant and R the Ricci scalar associated to the fixed worldsheet metric g"'^. The 
background charge Q is related to the Liouville parameter bby Q = 6 + 1/6, which provides the central charge 
cl = 1 + 6Q^. The matter part is described by a conformal action with a central charge cm = 1 — Qq^ where 
q = b — 1/6. The ghost part is the usual (6, c) system of the bosonic string theory, with central charge Cgh = —26. 
The vanishing of the total conformal anomaly leads to a constraint on the central chai^ges that fixes the Liouville 
parameter 6. In the following, we will concentrate on the p-critical models, realized as a LG with a matter sector 
given by a minimal model of the Lee- Yang series (2, 2p + 1) for which 6^^ = p + 1/2. 

For the critical gravity, the fields describing the conformal matter, metric and ghosts are formally decoupled. But 
in the coiTclation functions, any matter field needs to be dressed by the appropriate Liouville vertex operator to form 
a composite field of dimension (1, 1), $i je^'^^J'^^^^. The dressing charge Pij is related to the bare operator scaling 
dimension through 

where the indices (1, j) denotes the position of ^ij in the Kac table, and j = 1, • • • ,p. In order to respect the 
invariance under diffeomoiphisms, the LG composite fields should be integrated over the worldsheet. However, the 
presence of conformal Killing vectors allow to fix the position of some of these operators, provided we multiply them 
with a suitable ghost factor. To fix the notation, the integrated (or fixed) dressed operator with Kac label (l,p — J) 
will be denoted Vj with J running from zero to p — 1. 

On the disc, the Liouville action acquires a boundary term [13], 
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//B = -yZ/i cosh 7r6s (2.3) 



where K is the intrinsic curvature of the boundary.^ At the quantum level, the boundary states are labeled by the 
parameter s, which is related to the boundary cosmological constant ^b in the semi-classical limit by an hyperbolic 
cosine parameterization [14]. The Liouville gravity boundary states are tensor products of these states \s > with a 
Cardy state |(1, /) > describing the matter boundary conditions. Such states obey a linear decomposition property 

[15], 

i-i 

|s,(l,0>= Yl \s + iab, {1,1) > (2.4) 

a=~{l~l):2 

where the sum runs from — (/ — 1) to / — 1 with steps of two. This decomposition has been shown to be consistent 
with the matrix model boundary construction on the disc [12]. However, its validity for topologies of higher genus 
may still need to be checked more carefully [19,20]. 

Likewise in the bulk case, one may introduce dressed conformal matter boundary operators ^f .e^''-^'^^^^ with a 
dressing charge given by (2.2). The integrated dressed operators with Kac label (1, 2{l — k) — 1) where k is running 
from zero to / — 1 will be denoted as Bk- For convenience, one may consider the dressed operators as a perturbation 
of the LG action by introducing external sources, 

p-2 e-2 

AS = YtjVj + Y^ r'i^ Bk . (2.5) 

J=0 fc=0 



^In order to simplify the relation between the LG and MM boundary cosmological constants, the factor l/y^sm(-Kb^) has been put outside 
the definition of fiB, contrary to the notation used in [13]. 



The symmetry under translation of the Liouville field leads to associate a gravitational dimension to the coupling 
constants. Assuming a dimension two for the bulk cosmological constant /i, the other coupling dimensions are found 
to be [tj] = p+1 — J and [rf,] = I — k. It is noted that by definition, the bulk and boundary cosmological constants 

couple to the dressed identity operators, Tp_i = /i and r^_^ = fis, [fJ-B] = 1- The sphere and disc partition functions 
respectively scale as 2Q/b and Q/b. 

The bulk-boundary correlation number we are considering can be written as 

{VjB,)ll = ^(^+'=-'+V2)/2 ^ii) j^^p^^ ^^^ ,) (2.6) 

with Pj = Pi,p~j = ^(J + 1/2)6 and (3^ = l3i^2{i-k)-i = (^ ~ f^)b- Because of the presence of conformal 
Killing vectors, one may simply fix the position of both boundary and bulk operators, avoiding this integration over 
coordinates. The coefficient Afjl is independent of the Liouville parameter s and fi, it takes into account the matter 
and ghost factors, as well as some additional multiplicative constant in the expression of the Liouville correlator.^ 
This coefficient vanishes when the boundary operator Kac index does not belong to the fusion algebra of the modules 
associated to two copies of the bulk operator, i.e. 

M^ll = when I - k > p - J . (2.7) 

In the section 4.4 below, we will compare only the ^u- and s- dependent part of MM and LG correlators. In this 
manner, we do not need the explicit expression for the coefficients Afj^ and the property (2.7) will be sufficient for 
our purpose. 

The s-dependent factor in (2.6) has been obtained in Liouville field theory and is given in [17] as a Fourier 
transform 

fi(P.A»)=-,^d,e-»fi(P.A.), fl(P.A.) = ||±±|l|M p_3, 

where S± = a — ^/3 ^ P and S{x) denotes the double sine function [13] which satisfies the shift properties 

S{x + b) = 2sm{7rbx)S{x), S{x + l/b) = 2sm{-Kx/b)S{x), S{x)S{Q - x) = 1. (2.9) 

In the context of Liouville theory, the momenta P related to the charge of vertex operator by/3 = Q/2 — Pare 
purely imaginary. In this case, the Fourier transform in [17] is well defined and the contour C in the integral (2.8) is 
along the imaginary axis. However, in the context of Liouville gravity, these momenta are real and the convergence 
of the Fourier transform is no longer guaranteed. As a consequence, the contour C should be appropriately deformed 
(see details in the appendix A.l). 

One may notice that R{P, /3 , s) is an even function of P and s; 

R{P,P,s) = R{-P,(3,s), R{P,(3,s) = R{P,f],-s). (2.10) 

For the degenerate values j3 = kb,k ^ Z+, R is reduced to a ratio of sines due to the shift properties satisfied by S, 

RiPj, kb, a) = 4^-^ ntl sin(vrb^+ + .ab^) MnbS + vra6^) ^^.11) 

sin(7ro+/o) sin(7r()_/o) 

The contour integration can be performed, leading to an infinite sum over residues, 

R{Pj,kb,s) = . ,^ „ ,,, VTTsin{7r62(a + n)}sin{7r62(a + n) + 27rPj6}e 
sm{27TPj/b)\^^^J-^ 



-nbs{2n+k+2Pj/b) 



oo fc— 1 



(2.12) 



y 17 sm{7rb'^{a + n)} sm{TTb'^{a + n) - 27rPj6}e-^^"(2n+A,-2Pj/b) 



n=0 a=l 



^The exact expression of the factor Afjl depends on the proper normahzation of the operators. This difference of normalization between 
MM and LG correlators is sometimes referred as "leg factors", they coiTesponds to a multiplicative degree of freedom for the coupling 
constants. Here, it will be fixed by requiring monic polynomials in the resonance transformations of the section three. 



In the simple case k = I, the summation can be easily done and we get the simple expression 

jl(p b s) = 6 sinh27rPjg 

sm}i{2Tr Pj/b) sinhTrbs 

Unfortunately, the formula (2.12) is not convenient for comparison with the matrix model results. This is why in 
the next section we shall exploit the recursion relations over s satisfied by R in order to propose an alternative 
expression. 

2.2 Shift identities 

The bulk-boundary correlation R{Pj, kb, s) obeys several shift relations that are supposed to be related to the bulk 
[21,22] and boundary [23,24] ground ring structures. They provide a nice recursion relation which bypass the 
complicated summed expression and will turned to be useful below to relate MM and LG correlators. The readers 
who are not interested in the technical details might skip this subsection. 

The simplest relation is obtained from the expression (2. 12), when s is shifted by ±i/b; 

R{Pj, kb, s + i/b) + R{Pj, kb, s - i/b) = {- 1)'' 2 cos{2ir Pj/b) R{Pj, kb, s). (2.14) 

Another relation can be derived using the properties of the double sine function to show that 

R{Pj + 6/2, {k + 1)6, a) = 2{cos{2TTba) - cos 2Trb'^{k + Pj/b)) R{Pj, kb, a). (2.15) 

Plugging this into the Fourier transform, the first cosine just gives a sum over shifts of s, and we end up with 

R{Pj ± 6/2, {k + 1)6, s) = R{Pj, kb, s + ib) + R{Pj, kb, s - ib) - 2 cos 2Trb'^{k ± Pj/b) R{Pj, kb, s) (2.16) 

where a similar- relation was also derived for negative shifts. 

The most interesting relation is obtained if by shifting s and k but not Pj. The shift relation is properly derived 
in the appendix (A.2) and it reads 

cosh Tibs R{Pj, {k + 1)6, s) = 2 cos(27r6Pj) R{Pj, kb, s) - Ck [RiPj, kb, s + ib) + R{Pj, kb, s - ib)] 
smhTrbsR{Pj,{k + l)b,s)=isk [R{Pj,kb,s + ib) - R{Pj,kb,s - ib)] , (2.17) 

or equivalently 

cos2PjT:b „, „ , 1 /cosh7r6s sinh7r6s\ „,„ ,, 
R{Pj, kb, s±ib) = ^^R{Pj, kb, s)--[ ± i R{Pj, (k + 1)6, s), (2.18) 

Ck 2, \ Ck Sk J 

where we used the shortcut notation 

Cfc = cos7r6 k, Sk = sin7r6 k. (2.19) 

Together with the knowledge of the expression (2.13) for the operator Bi, this last relation fully determines R at any 
Pk- Indeed, it is easy to show that the following expression satisfies recursively the second shift relation in (2.17) 
with the initial condition given by (2.13) (see appendix A.3), 

^ '^' ^7=1 a=-(fc-l):2 /3=-(fc-i):2 

with Set = s + "ioib and Xq = x{sa) = ucosh-Kbsa- The expression for /c = 2 can be simplified into 

K{Pj,2b,s) - — . ,,„ p ,,s -2-— 2 1^^ ■ 1 , — • (^•^^) 

sm\i[2TTPj/b) cosh rrbs — c( ^ Sr±\ smh7r6s 



Finally, a recursion relation on k can be derived by considering R{Pj, kb, s + ib — ib), 
^Ck+iSk / 2 _ ^2^ j^^pj^ j^fj^ ^) ^ _^^ cosh7r6s R{Pj, {k + 1)6, s) + ^^^^°^ ^ ^ ~ ^^+^^ i?(Pj, (fc + 2)6, s) 

S2k+1 2s2fc+lSfc+l 



(2.22) 



with r = 2Pj/b in Cr- 



3 On the matrix side 

The one matrix model partition function is given by an integral^ over an A^ x A^ hermitian matrix M, 

^1MM= /dMe-f^^^(^) (3.1) 



where V{M) is a polynomial potential chosen to achieve the {p + l)-th Kazakov multi-critical point [25] (for a 
review of this model, see [1,2]). The partition function and the correlators can be expanded in powers of N^'^, each 
term being associated to a different topology. We focus here on the first term of the series which describe the planar 
topologies (sphere or disc). These first order terms depend on the t'Hooft parameter k^ = 1/g which weight the 
number of vertices of the planar Feynman diagrams. Since the parameter k controls the area of the discrete surfaces 
dual to the Feynman diagrams, it is sometimes referred as the bare cosmological constant. In the continuum limit, k 
is sent to a critical value k* where the mean area diverge, and we define the renormalized coupling e^'^ to = k — k* 
where e ~ iV~^/^P+^ is the lattice size cut-off. To this coupling is assigned the dimension [to] = p + 1. 

In order to compare the (p + l)-th multi-critical MM with the (2, 2p + 1) minimal Liouville gravity, we first 
need to introduce the KdV deformations [26-28]. These linear deformations of the potential V{M) by the other 

multi-critical potentials V{M) — )• V{M) + J2j \j ^v ^^^^ ^'^ ^^^ string equation that determines the string 
susceptibility - or two-punctured sphere partition function - u as a function of the deformations. Under a suitable 
normalization, the string equation reads 

uP+^-Y,tT^^^' = 0. (3.2) 

j=o 

where tj is the renormalized tj . The {p + l)-th multi-critical point, denoted *, is defined as the point 
where all the perturbations, including the MM cosmological constant to = ^o , are turned off, except for the 
coupling t -^ of dimension two. At this point, the MM can be compared to LG and t ^^ corresponds to the LG 
cosmological constant fi.^ In particular, n|* = ^/JI and the resolvent defined below is identified with the LG identity 
boundary Ipt function on the disc. 

As mentioned in the introduction, the bulk resonance transformation has already been treated in [7]. At the level 
of insertion of a single bulk operator, this unnecessary complication can be avoided by a linear redefinition of the 
polynomials that perturb the matrix potential. These redefined potentials, denoted here Vj, absorb the first order (but 
all orders in fi) of the resonance transformation. After this redefinition, the new potentials couple to bare parameters 



''This integral should be understood here as a formal series obtained by expanding the exponential, keeping only the quadratic term, the 
coefficient of which must be negative. 

Contrary to the unitary models (e.g. the 0{n) matrix model [29]), here the area is no longer measured by the dressed identity operators, 
but by the operator of highest dimension Vo- In this sense, the MM renormalized cosmological constant ij, does not corresponds to the LG 
cosmological constant tp_^ = /i. To avoid confusion, in the following we should only refer to fi under the denomination "bulk cosmological 
constant". 



tj that can be directly identified with the parameters tj perturbing the minimal LG in the continuum limit. In this 
setting, the bulk 1-pt functions are obtained from the matrix model correlators as, 



(tr Fj(M))cL = —g^ 



(K7)lg' (3-3) 



where the bulk couplings are turned off at *, except for tp_i = /i. With this choice, all the MM bulk 1-pt functions 
vanishes apart from the one with J = p — 1. 

Let us turn to the boundary effect. The disc with one marked point and trivial boundary conditions (i.e. leading 
to (1,1) matter BC in the continuum limit) is given by the resolvent, defined as 

W{x) = lim 1 (tr ^^\ = W'{x)+u;{x). (3.4) 

Af->oo A* \ X — M / 2 

This quantity is the Stielges transform of the eigenvalue density for the matrix M. The bare boundary cosmological 
constant x controls the length of the boundary of the discretized surfaces. In the continuum limit, it is also sent to a 
critical value x* where the mean boundary length diverges. This critical value can be taken to be zero by a shift of 
the matrix, and we define the renormalized boundary cosmological constant as ex = x, [x] = 1. In the process, we 
have to throw away the non-universal contributions which ai^e polynomials in x such as the term V'{x) appearing 
in the expression (3.4) of the resolvent. Since in the following we focus on the expression of the correlators in this 
continuum limit, we automatically throw away all the non-critical terms and write directly [30] : 

uj(x) = (tr ) =u' coshi-Ks/b), x = ucosh.TTbs. (3.5) 

\ ^-A^/cL 

The uniformizing parameterization of x has been introduced in order to resolve the branch cut over ] — oo, — n] of 
the resolvent. The right bound of the branch cut is identified with the "string susceptibility" u. 

On the disc topology, the boundary describing a matter with (1, Z) BC have been constructed in [1 1] by allowing 
the matrix to additionally interact with Gaussian vectors with flavors. For instance, the one and two-point functions 
of minimal LG with matter BC (1, 1) and cosmological constants x{s) (two identical boundaries for the two-point 
function) have been identified in [16] as 

where the leg factors have been included in the definition of the LG correlators and 

Fi{x,M)= TT {M — Xa), Xa = ucoslnrbsa, Sa = s + iab (3.7) 

a=-{l-l):2 

creates the corresponding MM boundary. The polynomials P^'->'^{x, M), determined in [16], encode the pure bound- 
ary resonance transformation and will be defined more precisely in the section 4.1 below. 

The MM correlators with identical boundaries Fi (x, M) can be derived from the following disc partition function 
in the presence of boundary sources ijj M", 

^«(i«) = /tr log (fi{x, M) + J] t(OMM \ . (3.8) 

To investigate the bulk-boundary matrix con^elator, we need to consider the disc partition function of Eq. (3.8) in the 
presence of the bulk sources in the CFT frame, as in (3.3). The matrix correlator we are interested in is defined as 






'"" dtjdt^' 



tr-——-trVjiM)) (3.9) 

Fi{x,M) / CL 



where 'c' denotes the connected part and * all the non-trivial couplings turned off. This quantity is easily evaluated 
from the expansion of the product Fi (see identity (A.l) of [12]), 

Ofaix) = - Yl n {^a-XprUt,——-tTVj{M)) . (3.10) 

We recursively make use of the property 

i-i /-I 

Y, n (2^a - a;/3)-'x^ = 0, fora<Z-l, (3.11) 

a=-{l-l):2 /3=-(i-i):2 

in order to rewrite the matrix correlators as 

i-i i-i / 1 \ 

0'jli^) = - E n i^a - xpr'x'^J tr -— j^tv Vj{M)) ^^ (3.12) 



a=-(l-l):2 /3=-(i-i):2 



c,CL 



This decomposition involves the derivative of the known bulk 1-pt functions with respect to the boundary cosmolog- 
ical constant [10], 

0^fi(.) = -(tr^trVAM)) = i,-V-V.-.5H||^ (3.13, 



X — M / c CL ^ ^^^'^ '^^^ 

where the normalization has been fixed such that 

0^^\ o(-) = -9.Ufixe.c. = l6-V-3/^ "^^-f/^-^) . (3.14) 

'^ 2 smnvros 

We end up with the following expression for the bulk-boundary matrix correlators, 

a=-(l-l) ;3=-(i-i) ^ "^ 

It should be noted that the s-dependent part of Ojq{x) exactly reproduces the correlator of LG, R{Pj, kb, s) in 
Eq. (2.20). 

Finally, let us stress that the expression (3.15) satisfies shift relations similar to those obtained for the LG coitc- 
latorin section 1.2. 

Ofai^ + i/b) + 0[;iis - t/b) = 2i-iy-'^coshi27rPj/b) 0«(.), 

'^0%+ii^) = 0%,^{s) + 0f_,^{s), (3.16) 

0f^{s±^b) = 0^j'~^%{s)-x^,0^;:'\s). 

The third identity is derived from a relation that can be linked to the insertion of a boundary ground ring operator, as 
explained in the appendix of [16]. 



4 Comparison of matrix models with Liouville Gravity 

4.1 The resonance transformation 

In order to compare the matrix model with the LG results, we have to take into account a possible redefinition 
of the coupling constants known as the resonance transformation. This transformation is a consequence of the 
ambiguity in the definition of the so-called contact terms, i.e. the values of the correlators that contain operators 
taken at coinciding points. Such contact terms can be reabsorbed into a finite renormalization of the couplings. 
Coinciding bulk operators lead to pure bulk coupling resonance. This phenomenon was first introduced in the 
matrix model context in [6] and then later investigated for the perturbed sphere partition function in [7]. When 
dealing with a worldsheet having a boundary, we have to take into account boundary operators at coinciding points 
which leads to pure boundary couplings resonance, i.e. redefinition of the boundary couplings involving only other 
boundary couplings. In addition, we also have to take into account a bulk operator coinciding with the boundary (or 
equivalently with its mirror image) and bulk operators coinciding with boundary ones. These two phenomena lead to 
a bulk-boundary resonance which translates into the presence of bulk couplings in the resonance transformation of 
the boundary ones. It is stressed that no boundary couplings can be involved in the resonance of the bulk ones. The 
resonance transformation can be seen as a finite renormalization and this should be related to the study made in [31], 
where it was noticed that no boundary couplings arise in the RG equations of bulk couplings. On the contrary bulk 
couplings in the RG equations of boundary ones generate induced boundary flows. 

In the problem we are considering, we have already taken into account the bulk resonance transformation in 
a suitable definition of the matrix potential deformations Vj, and the general boundary resonance transformation 
writes 

4' = ri') + E c,, n -j'rP '\ [4'^] = E PArj] + E ^'^^'^ (4.1) 

p,i^ J,k J k 

where the Cp j^ are just numerical constants. The resonance condition on the coupling constants dimension that 
restricts the form of the RHS gave its name to the transformation. We recall the dimensions of the couplings: 
[tj] = p + 1 - J and [r^/^] = I - k. 

The bulk and boundary cosmological constants fi = rp_i and fis = vJi are not turned off at the point * 
where the other perturbations vanishes. These two couplings are of minimal dimension among respectively the 
bulk and boundary ones, so that no other couplings can appear in their resonance transformation: ti_i = /i and 
tj_-i^ = ——fiB with Sfc = simrb'^k (or x = ^b, the coefficient —si/si appearing in the normalization of t^^-^ is 
due to the fact that the monomial of degree Z — 1 in Fi{x, M) is proportional to — f^x, see [11, 12]). The boundary 
resonance transformation relevant for the study of boundary one and two points functions, as well as bulk-boundary 
correlators, writes 

k=a k,m=a J=p+l—l+a 

The function P^' and Q'^' are polynomials in fi and fiB and the summations are restricted because of the resonance 
condition on the gravitational dimensions. Multiplication of the boundary couplings by a constant will only change 
the normalization of the boundary operator, leading to a different leg factor. In order to fix this degree of freedom, 
we assume here that at the first order ta = Ta , or P^ = 1 . 

To understand how the transformation (4.2) relates the correlators of the two theories, we first consider the LG 
boundary 1-pt function with a non-trivial operator inserted (B^ ^ Bi_i) and apply the chain rule. 



OTf^ a=o OTu Ota 



E-i-(-^L- 



a=0 ^ -'V--;/CL 



since the matrix correlators with powers of M", o < / — 1 are vanishing due to (3.11). In a similar way for the 
identity operator we get 



(^^-i) 



9Z(0 



(0 



a=0 '-'^l-l '-'^a. 



Y.d,sPil^ (tr ^)cL ^ (-l)'^^^)^'^-^^ (4.4) 



a=0 



,(0 



This is in indeed what we observe [12] since Pi'^^{fi, ^b) = —{si/si)^b- 

We now turn to the boundary 2pt functions and start again with the non-trivial operators A; / / — 1 / m, 



(BkBm) i^Q 




a=0 ^^m Ea , 



dr^^drli} d^ 



+ 2^2^ r. (l) ^ (l) ^(l) r.M) 

a=0 b=0 OTm OTf^ Ota Ot^ 



(4.5) 



dhP 



tr 



W 



dr?dTli}\ Fli^^M) 



CL 



EE 



ati" Btf 



a=0 fe=0 "^'m C'^fc 



tr 



^a+b 



W^JOX FKx,M) 



CL 



As in the case of the boundary Ipt functions, the first sum vanishes so that the quadratic term in the resonance 
transformation (4.1) plays no role here. The second term defines the two polynomials that were determined in [16], 

k 



{BkB,n)i^Q 



tr 



p(0,fc(a;^M)pW''"(x,M) 
Fi{x,My 



,(0,fc 



, P«''=(x,M) = ^Pi')''=(^,x)M". (4.6) 

I -I), 



CL a=0 

Since we imposed P^^'''''^ = 1, pW^^ is a monic polynomial in M. When one of the operator is trivial (m 
the computation is similar to (4.5) but we have to define the matrix polynomial separatly as 

i-i 



pm~i(^x,M) 



Sl 



Y,d.PP{fi,x)M\ 



(4.7) 



a=0 



where the normalization has been fixed in order to have again a monic polynomial in M. In the following we will 
not specify anymore whether or not the operators ai^e different from the identity, but just keep in mind this subtlety 
in the definition of p(0.'-i 

Finally, let us investigate the case of interest for this paper, the bulk-boundary correlators 



<''^-'>S - ^ 



,a=0 (^T^k '^'-a - 



h 5rf drjdtP 



+ E 



dH^ a^« 



.{I) «,(') 



(4.8) 



,=0 drjdr]^' dt'a> 

As in the case of the boundaiy Ipt and 2pt functions, the second term vanishes because either the correlator is zero 
(a < Z — 1) or a = / — 1 and the coupling tj_^ has no resonance involving tj. The first term in (4.8) is a sum of two 
contributions since the variables ta may depends on tj. 



{VjBk)f^G 



A dtP 522« 



.(0 



a=o drl dtjdta 



m 



^kl + J_P~l g^il)g^il) ^2^(;) 



5rP drj Qt^ Qt^ 



(4.9) 



Introducing the previous polynomials pW''"' and the bulk-boundary encoding polynomials^ 



(4.10) 



6=0 



fe=0 



*Note that since the relative normalization of the bulk and boundary couplings have already been fixed, there is no reason to require that 
^"''■^ is a monic polynomial in M. 



we finally get 

The gravitational dimension of the polynomials can be deduced from this relation, 

[{VjBk)] = k + J-l + l/2, [Ofj = a + J-l + l/2 ^ [P(')'^] = k, [Q(')"^] = J+l-p-l. (4.12) 

The formula (4.11) implies that we have to consider two different cases. When the LG correlator vanishes, both 
terms in (4.11) contribute. It appears that imposing the vanishing of the combination of the two MM coiTclators 
provides a sufficient number of constraints to fully determine the bulk-boundary resonance term Q^')''^. It will be 
done in the section 4.3. 

On the other hand, when the LG coiTclator is non-zero, i.e. l — k<p—J,-we observe that the term which involves 
qWJ is not present in (4.1 1). Indeed, in this case the boundary operator B^ does not belong to the fusion outcome of 
the two copies of the bulk operators involved in the resonance of ^fc<;_Lj_„_i- Because of the orthogonality property 
of the boundary 2pt functions, there can be no contribution from a contact term aiising when a bulk operator meets 
the boundary. We will explicitly check this result by relating MM and LG correlators for several non-trivial cases 
(A; = 0, 1, 2, / - 2, / - 1) in the section 4.4 below. 

4.2 Preliminary checks 

As the first verification of the formula derived for both MM and LG con^elator, we check the agreement of both 
expressions when no resonance is involved. This is the case for the boundary operator Bq since the boundary 
coupling Tq only enters the resonance transformation of Vq in (4.2) and pW'*^ simply equals to one. For the non- 
vanishing case I < p — J where Q^^'''^ is not involved, the expression (3.15) of Ojq{s) reproduces (2.6) for the LG 
correlator (Vj-Bo)lg' 

0,?>(,) = .-«/^<fl(P„i.,,), <' = -i,-32-- 1°™''' . (4.13) 

As the second verification of the previous formalism, we may investigate the shift relations involving s — )• s±i/b. 
More precisely, we should check that the sign flips arising due to such shifts (see (2.14) and (3.16)) are compatible 
with the formula (4.11) and the definition of the polynomials p(0.'«_ in particular, the coefficients Pa (/i, x) are 
polynomials of degree x'^"'^ in x. Since fi has a gravitational dimension two and x one, their monomials are of the 
form x'"""""^" where n is an integer. This implies Pa {fi, —x) = (— l)'^~"Pa {fJ-,x). Focusing again on the 
non- vanishing LG coiTclators, we have 



u 



j+k^l+y2j^() ^^^p^^ ^^ _ ^^^^ ^ ^ .^^^ ^ p^p^^ ^^ _ ^^^^ ^ _ .^^^^ 



' Jk 
a=0 K^-^'^) 



-{-\rH cosi2nPj/b) Y^ Pi')'^XM, x)0f^{s) 



as required. 
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4.3 Vanishing correlators and their constraints 

We consider the case I — k > p — J for which the LHS of (4.1 1) is equal to zero. It is easy to see by recursion over 
k that all the monomial terms have to vanish independently, resulting in / + J — p constraints 



tr 



Qm'J{x,M)M''' 
Fi{x,MY 



0]jI{x) for 0<a<l + J-p. 



(4.15) 



CL 



Solving these constraints allows to determine the polynomials Q^'■>'^^ {x,M) of degree l + J — p — linM that 
encodes the bulk-boundary resonance transformation. The bulk-boundary MM correlator Oj^ has been computed 
in the section 3. It is convenient to introduce the Chebyshev polynomials of the first and second kind, 



Tk{x) = u cosh/c7r6s, Uk{x) = u 
and to rewrite the expression (3.15) as 



f.smh.{k + l)7r6s 
sinh Tibs 



(4.16) 



ofs-) = {-ir'\h-^u^'~^^~'/' 



i~i 



i~i 



sm 



H^s/b) Y. n (^ 



-l„a Tp-j[Xa) 



a=-{l-l):2 ;3=-(!-l):2 



'sinh(7r6sQ,) 



l-i 



l-i 



ucos 



h{7rs/b) ^ Y\_ (Xa-X/s) ^x'^Up-j-i{Xa) 



a=-a-l):2 /3=-(i-i):2 



(4.17) 



where we used the expression of 2Pj = 1/6 — (p — J)b to expand the hyperbolic sine in the numerators. When 
a < I + J — p, the second sum in (4.17) vanishes because of the property (3.11) and we simply have. 



,(0, 



i' 



i-i 



i-i 



OZi^) = (-1)^-1-6-2^2^-^-1/2 sinh(vrV6) ^ H (^" " ^/^) 

a=-{l-l):2 l3=-(l-l):2 



-1 -^a P—Jv^o^l 

sinh(7r6sQ,) 



(4.18) 



The general boundary 2pt functions have been studied in [16], but a simplification occurs when the two bound- 
aries are identical. It can be seen by looking at the following correlator. 



tr 



M°- 



{M-y)Fi{x,M) 



CL 



E n (..-.«)-'< "<'''-'^<""' 



a=~(l-~l):2 /3=-(!-i):2 



y-xc 



(4.19) 



When y belongs to the set of {xa}a, all the terms vanish apart from Xa = V- For this remaining term, both numerator 
and denominator cancels and we end up with the resolvent derivative at Xa, 



tr 



M" 



(M-x„)FKx,M) 



i-i 



CL 



-6-V-V2 J] (.. - .,)-i.« ^^"^("^"/^) 



;3 = -(i-l):2 



sinh Tibsc 



(4.20) 



This result can be plugged in the expression of the MM correlators with two boundaries, to get 



tr 



Fi{x,MY 



l-i 



l-i 



CL 



(_l)'6-2^P-i/2 sinh(7r5/6) ^ J[ {x^- xp)' 



2 Xf^ 



a+b 



a=-(l-l):2 /3=-(l-iy.2 



sinh Tibsr 



(4.21) 
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Comparing (4.18) and (4.21), we propose an ansatz for the values of the polynomials Q^'^'"^ at the points M 



Xct- 

l-l 



Q«'-^(x,x„) = -In2(j-P) J] (x, - xp) [rp„j(x,) + g«'^(x,)sinh7r6s, 



(4.22) 



3=-(!-l):2 



where q^^^''^{xa) is a polynomial of degree p — J — linxa whose coefficients depend on s. This term arises because 
of the freedom in the identification to add x^sinhirbsa with c < / — 1 — a to Tp-j{xa) within the sum over a. 
Thus, the expression of (5'^'^'"^(x, Xa) contains p — J free parameters that will be determined a posteriori. Knowing 
/ values of Q^'^''^, we can use the Lagrange interpolation formula to obtain Q(')''^(x, M) as a polynomial of degree 

l-l, 

a=~{l~l):2 " 

But we know from scaling arguments that (^^'^''^(x, M) must be of degree / — 1 — (p — J), so we have to impose the 
vanishing of the p— J highest degree terms. These are exactly the number of constraints needed to determine q^^'''^ . 
These constraints can be derived from the study of the asymptotic at M — )• oo, 

l-l 

^ \Tp_j{xa)+q^^^'\xo,)smh.'nhso, xJJ = 0, n = Q---p-l- J. (4.24) 

a=-(/-l):2 

To solve these constraints, we first note that we can replace x^ in the previous equations by a Chebyshev polynomial 
of the first kind with the same degree. Then we decompose gW''^ over a basis of Chebyshev polynomials of the 

second kind, 

p-i-J 

g«'^(x„)= Y. i^''is)uP-'-'Uk{xa) (4.25) 

fc=0 

which allows to simplify the constraints 

l-l / \ 

y cosh(p — J)'Kbs a cosh nirbs a + 2^% sinh(/c + l)7r6sQ coshnvr^So, = (4.26) 

a=-(l-l):2 \ k ) 

and eventually perform the o-summation using 

l-l l-l 

y sinh77r6sQ = ^^ sinh77r6s, > cosh 77r6so = ^- cosh 77r6s. (4.27) 

a=-{l-l):2 ^ a=-{l-l):2 '^ 

We end up with a linear system over qj, , 



p-i—J ,j- -. -, 

Y^ 7/ ~{l)J ~'n 7/ I sr^ S{k+l±n)l , s _ i V^ S(p-J±n)Z , 



(4.28) 



that can be solved by inverting the {p — J) x {p — J) matrix U{s)J This technique provides a unique expression 
for the polynomials (5^'^'"^(x, M) that solves the constraints (4.15). From this expression, it is not obvious that the 
coefficients of (^('^''^(x, M) ai^e polynomials in x and /i. This can be shown, provided that the determinant oiU 
cancels with the numerators of each coefficients. This is indeed what happened in the few cases we checked. 



''For a matter of clarity, here the vectors and matrix indices start from zero instead of the usual convention one. 
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As a crosscheck for the expression (4.23) of Q*-'^'"^, we consider the insertion of a bulk identity operator, J = 
p — 1. In this simple case, there is only one constraint that can be easily solved, 

sinh irbs 

But since the cosmological constant /i = tp_i is not turned off at the critical point * and is not resonant, Q^'-'^P"^ is 
just the derivative of the boundary matrix operator Fi{x, M) with respect to /i at fixed x, 



Q«'^-Hx,M)=^^'("'^') 



dfi 



X fixed 



d cosh Tibs d 

dfi 2TTbfj, sinh nbs ds 



Fi{x,M) 



Fi{x,M) ^ smhiTTb'^a ^^'^^^ 



2-v/7Isinh7r6s ^-^ M — x„ 

^^ a=-{l-l):2 

This result is in agreement with the expression (4.23) and (4.29). 

As another example, let us investigate the case J = p — 2 where the matrix U is simply 2x2. The inversion 
gives 

~(0,p-2 ^ QCi 1 ,(/),p_2 ^ 2x^ - (2 + C2)xu'^ u^^ ^^ ^^^ 

x^ — n^c^ sinhyr^s' 2{x'^—u'^c\) sinhvrfes 

Even in the case of / = 4, the expression for Q*-^^'^"^ is rather complicated and we will not give it here. However, 
we have been able to check that the factor x^ — v?'c\ appearing in the denominator of the two previous coefficients 
cancels so that each monomial of Q^^''P~'^{M) is indeed a polynomial in x and fi. We believe that the cancelation of 
deiU is a general feature such that the polynomial Q^^^'P~'\M) determined by this method always has the required 
form. 

4.4 Non- vanishing correlators 

We investigate here the relation (4.11) in the case of non-vanishing LG correlators. As already mentioned above, 
there is no bulk-boundary resonance involved. The non-resonant case fc = has already been treated as a preliminary 
check. In this section, we restrict ourselves to the cases A; = 1, 2, / — 2 and I — 1. The aim is not to derive a general 
proof but only to provide convincing arguments for the consistency of the approach presented in section 4. 1 and the 
pure boundary resonance expression found in [16]. Accordingly, we will concentrate on the ^- and s- dependent 
part of the LG and MM con^elators which are already non-trivial. In particular, we neglect the coefficient Mjl that 
contains the matter and ghost sectors contributions. With a more careful analysis, one could also a priori derive their 
expression from the MM side but this is beyond the scope of this paper. 

The pure boundary resonance transformation, encoded in the polynomials P^')'''"(x, M), has been determined in 
[16] by solving the orthogonality constraints on the boundary 2pt functions. Unfortunately, their explicit calculation 
is tedious when k is large since it involves the inversion of a matrix of size (/c + 1) x (/c + 1). This is why in appendix 
B we developed an alternative derivation which is useful for k close to /. 

Case k = 1: The resonance polynomial is of degree one and reads [16], 

P(')'^ (x, M) = M-^— (4.32) 

Q-i 

To compute the MM correlator involving p(')'^(x, M), we need the expression of O j[. The later can be obtained 
using the third equation of (3. 16) to express it in terms of shifted Ojq, 



i(0^.^ _ 1 (ni^~^)(. ^ ,h\ ^ roC-i)/. _ .•/,^^ ^ ... .rnW, 



<^jM = o [^V\^ + ib) + 0%'\s -ib)]+ xci.^O^^lis). (4.33) 
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The sum of the two shifted MM coiTelators can identified with LG coiTelators through (4.13), and simplified using 

(2.17), 

2 2 

0%^\s + ib) + 0[j-^\s - ib) = -A'-^-^^u-^~'+''/^Mf^R{Pj, (l - 1)6, s) + 2x^^O%) (4.34) 

^l—l Q — 1 

where we denoted r = 2Pj/b. We deduce 

(tr^^g^try,(M)^ =n^-^^^/^MflRiPMl-ms), AA« = -2^<). (4.35) 

Case k = 2: The expression of p(')'2 is known explicitly (see the formula (4.13) of [16]), 

P(')'2(x, M) = m2 - 2^xM + (£bi±5z^^2 _ fifk^2_ (436) 

Q-2 ■S2i-3Q-2 ■52^-3 

To avoid too complicated expressions, we simply show the proportionality of the MM correlators with R{Pj, {I — 
2)6, s) but the coefficient A/j2 could also be determined. We notice that the insertion of the polynomial 

P(^^'\x, M) = /'^^'^|. = M2 - 2q„ixM + x2 - uhU (4.37) 

within the MM correlator trivially gives an expression proportional to R{Pj, {I — 2)6, s). Consequently, we only 
need to show that the difference of the two polynomials, now of degree one in M, is also proportional to this LG 
correlator, 

P«'2(x,M) = xM - ^'i~2x' + {ssi., + siW ^ pm^^^M) - P«'2(x,M). (4.38) 

4s2Z-3 

It can be rewritten by introducing the polynomial p(')'^ of (4.32), 

P«'2(x,M) = xP«'i(x,M) + ^l^^fl^{x^-u^cU), (4.39) 

S21-3C1-1 

and inserted into to the MM correlator to give a sum of two LG correlators through (4. 13) and (4.35): 



F,{x,M) '1^ '"c-i 



-2xucrR{Pj, {I - 1)6, s) + —^ ^^^R{Pj, lb, s) 

S21-3S1-1 

(4.40) 

The expression inside the parenthesis can be simplified using the recursion relation (2.22), it is indeed proportional 

to R{Pj, {I - 2)6, s). This shows the vahdity of (4.1 1) with k = 2. 

Boundary identity operator (k = I — 1): All the bulk-boundary LG coiTelators containing the identity boundary 
operator -B/_i are non-vanishing, regardless of the bulk operator. The simplest way to determine the polynomial 
p(0,'-i is to use the fact that the boundary cosmological constant is not turned off at the *-critical point. Since this 
coupling is not resonant, p(v.'-i js simply the derivative of Fi with respect to t\_-^ = —f-x, 

P^^^^^-\x,M) - ^^ ^-^'(^'^) _ ^1 Fii^^M) y. sinhTTbSa ^^^^^ 

Si dx Si sinh nbs ^-^ M — x„ 
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Pluging this result into the expression (3.15) for the MM correlator and using (4.27) to perform the a-summation, 
we obtain 

(4.42) 
The expression of A/j/_]^ at Z = 2 is in agreement with A/j/ given in (4.35). 

Case k = I — 2: The polynomial p(0.'-2 ^^^ ^^ computed using the method given in appendix B. The matrix U 
to be inverted is only of size 2x2, and we find 

.m-2 _ 2S21XU ^(0.^-2 _ siu^ (443) 

" S2A(x) sinhyrfts' ^ siA(x) sinhyrfts' 

with the determinant 

A(x) = detZY = 4^M^^ (^2 _ ^2^2N ^ (444) 

C1S2S3 

Using these results, we can derive the MM correlator associated to the polynomial p(')''~^(x, M), 



p(^).'-^(. M) ^^ ^^^^^ \ ^ ,.-3/2^W ^^(p^^ 26, . 



tr — -— ^-^— ^tr Vj{M) = u^-'/'N'ji'RiPj, 26, s) (4.45) 



F,(x,M) -^^ Y 



c,CL 



with the expression (2.21) for the LG part and (r = 2Pj/b) 

AA,,_, = --6 Bmh(2.P./6) ,^,^,^(,2_,2)(,2_,2) • (4-46) 

This coefficient is in agreement with the formula (4.13) for / = 2 and (4.35) for / = 3. 

5 Summary and discussion 

In this paper we investigated the bulk-boundary disc correlators of the hermitian matrix model and the ( 2 , 2 p + 1 ) 
minimal Liouville gravity. In our study of LG, we specialized to the case of a coupling to degenerate matter operators 
for which the Liouville dressing chai^ges takes only a finite set of values. We derived two different expression, (2. 12) 
and (2.20), for the cosmological constants dependent part of the correlator. We also provided various useful relations 
involving either a shift of the boundary parameter ((2. 14),(2.16) and (2.17)) or a recursion relation on the boundary 
operator momentum (2.22). Given the large range of applications for the Liouville theory, such relations could also 
be useful in other contexts. 

In the second section, we constructed the MM bulk-boundary correlator following the method of [11, 12]. We 
obtained its expression in the continuum limit (3. 15) and successfully identified it to the LG correlator in the non- 
resonant case. But our main interest was in the study of the resonance transformation at the bulk-boundary level. 
For this purpose, we considered this transformation in a very general setting in section three. It led us to distinguish 
between two different cases. In the first case, the bulk-boundary LG correlator is vanishing due to the fusion rules 
obeyed by its matter part. Imposing this vanishing condition on the MM correlators, we were able to find the first 
order of a bulk-boundary resonance (4.22). This transformation can be encoded in a matrix polynomial whose 
coefficients satisfy a the Unear system of equations (4.28). 

In the second case, the LG correlator is non-vanishing. There, no bulk-boundary resonance arises - except for 
the identity operator - and the MM and LG correlators were seen to agree in a number of specific cases, provided 
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we take into account the boundary resonance investigated in [16]. In this process, we derived an alternative method 
to compute the matrix polynomials encoding the transformation which is particularly efficient for the boundary 
operators having a small Kac index. 

The explicit agreement between LG and MM correlators in the non-vanishing case has been shown only for some 
particular cases. Even though these cases are enough to draw the general pattern, a complete proof is still lacking. 
Such a proof would provide the expression of the factors Mjl that can be identified, once a proper normalization is 
fixed, to the matter contribution of the LG correlator, namely to the bulk-boundary disc correlation function of the 
(2, 2p + 1) minimal model. 

Our ultimate purpose is to conjecture the expression of the resonance transformation including the boundary 
effects at all order. Such a proposal was made for the bulk case in [7] after the study of the 3-pt functions. The 
boundary disc 3-pt function can be investigated in a similar way. In particular, the LG correlator is known and obeys 
the property of factorization into ghost, matter and Liouville sectors. The expression of the Liouville correlator 
found in [32] is rather complicated but should simplify in the case of a coupling to degenerate matter operators. 

We defined the MM boundary perturbations to be simple power of the matrix, in contrast with the bulk case 
where these perturbations are critical potentials related to the underlying KdV hierarchy. There may exist a better 
definition of the perturbations in the boundary case, possibly related to an underlying integrable hierarchy. This 
change of basis for the MM perturbation may allow to rewrite the resonance transformation in a simpler form. 

The relations we found in section 1 .2 for the LG correlators are believed to be consequences of the presence 
of the ground ring structures and the connection still need to be specified, as we did for the boundary 2-pt function 
in [16]. The exact reahzation of such structures within matrix models is still an open problem. For instance, in the 
0{n) model the boundary ground ring identities have been related to the continuum limit of MM loop equations 
in [33]. At the moment, no general picture including these two examples has emerged. 

Finally, we could also apply this boundary construction to other matrix models, such as the 0{n) [34], ADE [35] 
or dimers [36] models. These models have an interpretation as statistical models defined on a fluctuating lattice. We 
hope to develop a statistical interpretation of the matrix boundaries in this context. 
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A Technical details for the LG side 

A.l Fourier transform 

The asymptotic behavior of the LG correlator when s — )• oo can be deduced from a scaling argument^ 

R{Pj, kb, s) ~ e^^'^o", ao = 2Pj -kb€R. (A. 1) 

The Fourier transform of R{Pj,kb,s) properly exists only when ctq < 0, i.e. when 2Pj < bk, and we first 
consider this case. Then, the Fourier transform can be inverted by choosing any straight contour lying in the strip 



*The LG correlator obeys a scaling relation that allows to write it in the form R{Pj, kb, s) = fig" F{iJ,/fi%). In the limit /i — )> 0, the 
correlator must remain finite so that -F(O) is a non-zero constant. In the algebraic limit s — >■ cxd, or /ib » y/J^, F can be expanded around 
zero and we deduce the asymptotic behavior (A.l). 
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Re a G]cro,-cro[, 

/•c+ioo 

R{Pj,kb,s)= dae-^^'''R{Pj,kb,a), c£]ao,-ao[ (A.2) 

Jc—ioo 

and the Fourier transform R{Pj, kb, a) is analytic in this strip. We now examine the poles of the expression (2. 11) 
for R. Since some of the zeros in the denominator cancel with zeros in the numerator, we are left with two infinite 
series of poles centered on iPj, poles being separated by a distance b. In addition, the two intervals of length kb, 

]±Pj- kb/2, ±Pj + kb/2[ are free of poles, 

a = ±Pj + -kb + nb, n>Oorn<-k. (A3) 

When o-Q < 0, the two intervals overlap and the strip Re cj g] — cjo, o-o[ around the origin is indeed free of poles. 
Closing the contour on the right, the sum of residues leads to the formula (2. 12) for the LG correlator. When (Tq > 0, 
the two intervals free of poles does not overlap anymore. However, we can still make sense of i? as a Fourier 
transform by smoothly deforming the contour from the case ctq > in the inversion formula. This contour, denoted 
C, is such that only the poles of (A.3) with n > are picked up. The residue formula again lead us to the expression 
(2.12) for the LG correlator. We can check that this expression has the correct asymptotic (A.l) and satisfies the two 
reflection properties, 

R{Pj, kb, s) = Ri-Pj, kb, s), R{Pj, kb, s) = R{Pj, kb, -s). (A.4) 

A.2 Derivation of the shift equation 

To derive the two shift relations (2.17), we first notice that 5± is invariant under k^k + l, a^a + 6/2 and since 
the shift of k increment the number of terms in the numerator product of (2.11), we get 

R{Pj, {k + 1)6, a) = 4sin7r6(CT + {k - 1)6/2 + Pj)sin7r6(o- + {k - 1)6/2 - Pj) R{Pj,kb,a- 6/2) (A.5) 

This relation can be inserted within the Fourier transform, and the shift over a can be reabsorbed under a change of 
variable a ^ a — 6/2. This change of variable shifts the contour of 6/2 to the left, so that we could a priori pick up 
an extra pole. However, since the poles of -R in o" are separated by a distance 6 and J, k are integer, the contour C can 
always be chosen such that no additional poles are picked up. This choice will be confirmed below by the reflection 
properties Pj ^ —Pj of the recursion relation, and we can safely write 

RiPj, {k + l)b, s) = -ie-^^' I dae''^''"' 4sin7r6((j + A;6/2 + Pj) simrb{a + kb/2 - Pj) R{Pj, kb, a). (A.6) 

Then, we use a trigonometric identity to transform the product of sines into a difference of cosines, 

R{Pj, {k + 1)6, s) = ie~^^' I dae"^^'"' Uinba^inb^k ^ ^~2tnba^~inb^k _ 2cos27r6Pj] R{Pj,kb,a). (A.7) 

Inside the brackets, the first two terms can be absorbed by a shift of the variable s, the third one does not depend on 
the integration variable, so we end up with 



2l.„/„ ,, ., N _i^b2^ 



e'^'^'RiPj, {k + 1)6, s)=2 cos 2^6Pj R{Pj, kb, s) - e'^" "RiPj, kb, s - ib) - e''^" "RiPj, kb, s + ib) (A.8) 

Note that this relation is invariant under the bulk reflection Pj ^ —Pj which confirms the fact that no extra poles 
were picked up since individual residues does not obey this symmetry. Finally, using the property R{P, /3, —s) = 
R{P, p, s) we easily derive the two equations (2.17). 



17 



A.3 Solution of the second shift equation (2.17) 

To show that the expression (2.20) satisfies the recursion relation given by the second equation of (2.17), we just 
plug it in the RHS of this equation and then shift the indices a to get 

D/D n , ^\■i. \ fc-i rfcsinhJTrfe^fc .^ ^^ -pr ^ sinh27rPjga±i 

± a=_(fc-l):2 /3=-(fc-i):2 "±^ 

A._i rfoSinhivrft^A; / ^-^ -p-i- i sinh27rPjSQ, 



smhyrfts \ -^^ -'■-'- smhyrfts 

\Q=-(fc-2):2 /3=-(fc-2):2 

E~ -TT . x-i sinh27rPjSQ 



a 



sinh T^hsr 
I 



a=-k:2 /3=-fc:2 " 



where we denoted 

2'=-^6 



fc-ii ^=^-1 



JJ sinh^ ivrfe^^. (A. 10) 



sinh(27rPj/6) J-i 

7=1 

We first examine the term a = k appearing in the first sum and see that in order to form a complete product over /3 
running form —k to k with steps of two (excepting a), we need an extra factor {xk — a^-A,)- Similarly the same extra 
factor is needed for the term a = —k. Then, all the other terms have contributions of both sums which differences 
leads to the same factor: 

H [Xa-Xp)'^- JJ {Xa- XpY'^ = {Xk-X_k) H {Xa-Xp)~'^ (A.ll) 

,3 = -(fc-2):2 ^ = -fe:2 /3 = -fe:2 

Thus, the correlator with shifted /? is equal to 

i?(p,,(fc+i)6,.)= ''''". (xk-x.,) Y. n (^"-^/^)"' ''lhlf" 



a=-fc:2 /3=-fc:2 



_i sinh 27r PjSq, 
sinh irbsa 

a=-k:2 P=-k:2 ° 



(A. 12) 



^'Tfe+i ^ n (2;q - 3;/3) 



B Alternative solution for the pure boundary resonance encoding polynomial 

The solution found in [16] for pW'*^ is manageable when k is small since it involves the computation of determinants 
of size (fc + 1) X (A; + 1). The alternative solution we demonstrate here is convenient when I — k is small, relying 
on the inversion of a matrix of size {I — k) x {I — k). However this solution is specific to the case of two identical 
boundaries. The orthogonality conditions arising from the boundary 2pt function can be written as [16], 
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Together with the constraint that P^^>'^{x^ M) is a monic polynomial of degree k in M, these requirements fully 
determine pW''^^ The MM correlators appearing in (B.l) can be computed using equ. (4.21), 



E n (^--^/^r^ JwL"^" =Q- y^<a<k. (B.2) 



a=_(i_l):2 /3=-(i-i):2 " 



The following ansatz for the values of P^^^'''{x, Xa) solves the orthogonality constraints, 



i-i 



P(')'^(x,x«) = Yl {xa-x,3)smh-KbsaP^^'^'''{xa), (B.3) 

I3 = -[1-1):2 

where p^'->'^{xa) is a polynomial of degree I — 1 — k in Xa with coefficients depending on x to be determined. Using 
the Lagrange interpolation formula, we define the polynomials p(')''^(x, M) as 

P(0,fc(,, M) = Fiix, M) ^ ^^^i^^^/)'^(x.). (B.4) 

These polynomials are of degree I — 1 but from scaling arguments pW''^ should be of degree k. Thus, we have to 
impose the vanishing of the l — k — 1 first coefficients. These l — k — 1 new constraints are not sufficient to determine 
the polynomials p^^^'^{xa) depending on / — fc coefficients. The remaining constraint comes from the condition 
that pW.'s is monic which also restricts its term of degree k. All these constraints can be read from the asymptotic 
behavior at M — )• oo, 

l-i 

y~^ x^ ^^'^'''(xq) sinhTT^So = 5n,i-k~i, 0<n<l - k. (B.5) 

a=-{l-l):2 

To determine p^'^'^, we first note that in the previous equation xJJ can be replaced by any monic polynomial of the 
same degree and we choose (suitably normalized) Chebyshev polynomials of the first kind. Then, we decompose 
p{i),k on a basis of Chebyshev polynomials of the second kind, 

i-k-i 

/)''(^a)= E P^J;}'" U-^'^''-'^"''^U^{Xa) (B.6) 

m=0 

in order to rewrite (B.5) as^ 

l-k-l l-l 

'^ Pm''^ X] coshmrbsa smh{m + l)7Tbsa = 2^~''~'^6n,i-k-i, 0<n<l-k. (B.7) 

m=0 a=-(l-l):2 

(I) k 

The a-summation can be performed using (4.27) leading to a linear system of equations of ^m that can be solved 
by inversion of the matrix U of size (I — k) x (I — k): 

l—l—k ^/_^._2 1 

V UnUx)p^J!'Hx) = -— --5„,;_fc„i, Un„^{x) = - V "^"^+^^"% ^±„(x). (B.8) 

^ smh irbs 2 ^ Sm+i±n 

m=U ± 



'Due to normalization issues for the Chebyshev polynomials, this expression is only valid for A: < Z — 1. When k = I — 1, the RHS is 
(5„,o instead of 2^^5,i,o- 
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This matrix U is the same that appeared in (4.28) for the determination of the bulk-boundary resonance encoding 
polynomials Q(')''^. 

As a crosscheck, we can consider the case of the insertion of the boundary identity operator k = I — 1. The 
relation (B.5) supplies only one constraint that can be solved, leading to 

pm-^ = ^— ^ (B.9) 

Si smnvros 

which indeed gives a polynomial p(0:'-i that corresponds to the expression (4.41) obtained as a derivative of Fi 
with respect to x. As another crosscheck, the explicit expression for the polynomial with / = 4 and k = 2 has also 
been derived, it matches the formula (4.13) given in [16]. 
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